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Abstract 

We address the problem posed by the inhomogeneous trapping fields when using ultracold fermions to simulate 
strongly correlated electrons. As a starting point, we calculate the density of states for a single atom. Using 
semiclassical arguments, we show that this can be made to evolve smoothly towards the desired limit by varying the 
curvature of the field profile. Implications for mutually interacting atoms in such potentials are briefly discussed. 
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Since the first observation of Bose- Einstein conden- 
sation in a trapped gas of ultra-cold atoms [1], there 
has been an explosion of experimental and theoreti- 
cal work on such systems. More recent developments 
include adding to the trapping fields a laser standing 
wave (a so-called 'optical lattice'). This generates an 
effective potential for the atoms which is given by 



/ 2nx\ 

V(x) = V COS [ — — I + Vtrap(z), 



(1) 



where a is the spatial period of the optical lattice, Vo its 
strength, and Vt Ia ,p(x) the effective potential due to the 
trapping fields. (The generalisation of (1) to more than 
one dimension is straightforward.) Recently, degener- 
ate fermions have been loaded into such optical lattices 
[2], opening up the possibility of simulating strongly 
correlated electrons. Here we explore this prospect by 
considering, as a starting point, the single-atom den- 
sity of states (DOS). 
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Fig. 1. Confining potential for various values of n, as indicated. 
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Perhaps the simplest case is that in which the optical 
lattice is very strong. In this case, we can replace the 
full Hamiltonian, 

V 2 ( 2ttx \ 

Hfull = 7^— + Vo COS + VtTB. V {x), (2) 

Am \ a J 

by a tight-binding version, 

H = ^2{-t\j)(j + l\-t\j + l) (j\ 

j 

+V^ P (ja)\j){j\}, (3) 

where \j) is a state in which the atom is in its local 
ground state on site j of the lattice. 

In models of electron systems, one usually assumes 
the confining potential to be a 'box': Kra P (a;) = for 
\x\ < L, oo otherwise (the solid line of Fig. 1). In the 
bulk limit, 

L -> oo, (4) 

Vtra,p(x) vanishes everywhere and one obtains the usual 
tight-binding DOS (the solid line of Fig. 2). 

On the other hand, ultracold atoms are usually sub- 
ject to a harmonic trap: Vtrap(x) = t(x/L) 2 . The bulk 
limit (4) corresponds again to a vanishingly weak trap- 
ping potential, yet the DOS is qualitatively different. 
It is given by the dotted line in Fig. 2, as can be shown 
either via a WKB approach [3] or a 'local-density ap- 
proximation' (LDA): 



P(e) = Phom(e - Krap(jo)), 



(5) 



where phom(e) is the homogeneous DOS shown in 
Fig. 1. (We prove elsewhere that these two methods 
are equivalent.) 
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This state of affairs is unsatisfactory, in the sense 
that one would like to be able to interpolate smoothly 
between the harmonic scenario realised in atom-traps 
and the open boundary conditions of usual solid-state 
physics — particularly if one hopes to use the former 
to simulate the latter. We propose to achieve this by 
varying the form of the power law at x = 0. Consider 
the family of potentials 



V tl!i p{x) 



(6) 



These reduce to the two cases studied above forro — > oo 
and n — 2, respectively. Let us therefore study the be- 
haviour of the single-atom DOS of (3), with trapping 
potential (6), as a function of n. We call this progres- 
sion through the family (6) 'finite-curvature scaling', 
by analogy with the finite-size scaling procedure used 
in numerical simulations [5] . 

We may obtain some information via straightforward 
asymptotics. Firstly, in the case e 2t, we may neglect 
the first two terms in (3) . What remains is just the T4ra P 
term, the eigenstates of which are clearly the position 
eigenstates \j). The 'dispersion relation' is then e = 
t\x/L\ n , implying that \x\ = L(e/t) 1 ^ n , from which the 
DOS may readily be obtained: 
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This gives e _1//2 behaviour in the harmonic case [3] in 
agreement with recent experiments [4]. Note that it 
also predicts e~ behaviour in the n — > oo limit; this 
prediction, though correct, is irrelevant, because the 
spectral weight at e > 2t goes to zero as n — > oo. 

Secondly, we may look at the case e w —2t: these are 
the states very near the bottom of the band. For them, 
the tight-binding form of the DOS can be expanded 
around k — 0: £hom(fc) = — 2tcos(ka) w — 2t + ta 2 k 2 , 
corresponding to the kinetic energy of a free particle 
with mass m e g = h 2 /2ta 2 . The total energy (classi- 
cally) would therefore be given by 



-2i + p 2 /2m cff +VtT& p (x); 



(8) 



using this relation as the orbit equation, and applying 
the WKB approximation as in the harmonic case [3], 
one can show that the DOS near the bottom of the 
band has the form 



p(e) ~ (e + 2t) 



(l/n)-(l/2) 



(9) 



This produces the correct square-root singularity in the 
n — ► oo limit, and the previously calculated constant 
behaviour in the harmonic trapping case. Notice that, 
according to (9), a singularity exists at the bottom of 
the band for any n > 2. 

To explore what happens at the top of the band (e ~ 
2t), we evaluate the DOS in the LDA approximation 



(5) numerically. The results for a few different values 
of n are shown in Fig. 2. 
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Fig. 2. The density of states for various values of n, as indicated. 
As n increases, the singularity at e — — 2t strengthens, while 
the one at e — It becomes increasingly asymmetric. 

In summary, we have presented an examination of 
the DOS of a single atom in a strong optical lattice plus 
a shallow power-law confining potential. Varying the 
power law permits one to scale from the case of har- 
monic trapping to the hard- wall case relevant to solid- 
state physics, and we call this variation finite-curvature 
scaling. The next step, for fairly weakly interacting 
atoms, would be to introduce the atom-atom interac- 
tions perturbatively. In the fermionic case, one can see 
that something special will happen (enhanced suscep- 
tibilities, or perhaps even ordering) when the chemical 
potential lies near the singularity at e = 2t. Analysis 
of this case is the subject of current research. 
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